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The singularity of the solutions obtained before in the teleparallel theory of gravitation is 
studied. Also the stability of these solutions is studied using the equations of geodesic devia- 
tion. The condition of stability is obtained. From this condition the stability of Schwarzschild 
solution can be obtained. 



1. Introduction 



Static uncharged black holes in general relativity are described by the well-known Schwarzschild 
solution. Singularity can be happen when a gravitational collapse takes place and continue 
until the surface of the star approaches the Schwarzschild radius, i.e., r = 2m [1]. Hawking 
and collaborators discovered that the laws of thermodynamics have an exact analogues in 
the properties of black holes [1]~ [4]. As a black hole emits particles, its mass and size 
steadily decrease. This makes it easier to tunnel out and so the emission will continue at 
an ever- increasing rate until eventually the black hole radiates itself out of existence. In the 
long run, every black hole in the universe will evaporates in this way. 

The tetrad theory of gravitation based on the geometry of absolute parallelism [5]~[14] 
can be considered as the closest alternative to general relativity, and it has a number of at- 
tractive features both from the geometrical and physical viewpoints. Absolute parallelism is 
naturally formulated by gauging space-time translations and underlain by the Weitzenbock 
geometry, which is characterized by the metricity condition and by the vanishing of the cur- 
vature tensor (constructed from the connection of the Weitzenbock geometry). Translations 
are closely related to the group of general coordinate transformations which underlies general 
relativity. Therefore, the energy-momentum tensor represents the matter source in the field 
equation for the gravitational field just like in general relativity. 
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The tetrad formulation of gravitation was considered by M0ller in connection with at- 
tempts to define the energy of gravitational field [15, 16]. For a satisfactory description of 
the total energy of an isolated system it is necessary that the energy-density of the gravita- 
tional field is given in terms of first- and/or second-order derivatives of the gravitational field 
variables. It is well-known that there exists no covariant, nontrivial expression constructed 
out of the metric tensor. However, covariant expressions that contain a quadratic form of 
first-order derivatives of the tetrad field are feasible. Thus it is legitimate to conjecture that 
the difficulties regarding the problem of defining the gravitational energy-momentum are 
related to the geometrical description of the gravitational field rather than are an intrinsic 
drawback of the theory [17, 18]. 

In an earlier paper [19], the author used a spherically symmetric tetrad constructed 
by Robertson [20] to derive three different spherically symmetric space-times for the cou- 
pled gravitational and electromagnetic fields with charged source in the tetrad theory of 
gravitation. One of these, contains an arbitrary function and generates the others. These 
space-times give the Reissner Nordstrom metric black hole. Calculations of the energy as- 
sociated with these black holes using the superpotential method given by M0ller [15] have 
been done [19]. It has been shown that unless the time-space components of the tetrad field 
go to zero faster than 1 / y/r at infinity, one got different results for the energy. 

It is the aim of the present paper to study the singularity of the three black hole solutions 
obtained before [19] and then derive the condition of stability using the geodesic deviation 
[21]. This study is important to gain more investigation about the solutions obtained before 
[19]. In §2, a brief review of the three black holes are given. The singularity problem of 
these black holes is studied in §3. In §4, the condition of stability is given. Final section is 
devoted to main results. 



2. Spherically symmetric black hole solutions 



In a previous paper the author used the teleparallel space-time in which the fundamental 
fields of gravitation are the parallel vector fields by} 1 . In the Weitzenbock space-time the 
fundamental field variables describing gravity are a quadruplet of parallel vector fields [22] 
b^, which we call the tetrad field in this paper, characterized by 

D u b^ = d u b^ + r M A jA A = 0, (1) 

where T^\ v define the nonsymmetric afline connection coefficients. The metric tensor g^ is 
given by 

9/iu b pbtvi 

where summation convention is taken over i. Equation (1) leads to the metricity condition 
and the identically vanishing curvature tensor. 

The gravitational Lagrangian L G is an invariant constructed from and the contorsion 
tensor 7^ given by 

1/j.up b pb; w -^ p — {Typp + Tpp U Tpvp) Tpvp b pbi V ^ p b yb%p^ p, (2) 



where the semicolon denotes covariant differentiation with respect to Christoffel symbols and 
comma is the ordinary differentiation and T pup is the torsion. It is of interest to note that 
the concept is as old as the gravitation theory of Einstein. The torsion notion of a variety, 
besides the curvature was introduced by Cartan [23] that also gave a geometric interpretation 
for both tensors. In teleparallel theories the gravitational interaction is described by a force 
similar to the Lorentz force equation of electrodynamics, with torsion playing the role of 
force [18]. 

The most general gravitational Lagrangian density invariant under parity operation is 
given by the form [7, 22, 24] 

C G = ^g~L G = + a 2l ^ p lpvp + a^ vp lpvp ) (3) 

with g = det(g pu ) and $ M being the basic vector field defined by $ M = l p pp - Here 0:1,012, 
and 0:3 are constants determined such that the theory coincides with general relativity in 
the weak fields [15, 22]: 

«! = --, a 2 = -, o 3 = -(l-A), (4) 

K K K 

where k is the Einstein constant and A is a free dimensionless parameter*. 
The electromagnetic Lagrangian density L e . m . is [25] 

Le.m. = - l -g w g v °F pv F pa , (5) 

with F pu being given by 1 " F pv = d p A u - d u A p . 

The gravitational and electromagnetic field equations for the system described by L G + 
L e m . are the following: 

G pv + H pu = -kT^, = 0, d v (V=gF^) = 0. (6) 

with G pu being the Einstein tensor of general relativity defined by 

G pu = R pu — -g pu R, (7) 
R^({}) is the Ricci tensor defined by 

R,M) = d P {%} - d v {^} + {Q {£,} - {U 

and R({}) is the Ricci scalar 

R({}) = g^R pu 

H^ v and K pu are defined by 



Hpv = A 



(8) 



*Throughout this paper we use the relativistic units, c = G = 1 and k = 8tt. 
^Heaviside-Lorentz rationalized units will be used throughout this paper 
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and they are symmetric and antisymmetric tensors, respectively. The energy-momentum 
tensor T pv is given by 

= -g pa F^F ua + ^ u F pa F pa (10) 

It can be shown [22] that in spherically symmetric case the antisymmetric part of 
the field equations (Eq.(6)) implies that the axial- vector part of the torsion tensor, a p = 
{l/2>)e lJlV p a r ) upcr , should be vanishing. Then in Eq. (8) vanishes, and the field equa- 
tions (Eq.(6)) reduce to the coupled teleparallel equivalent of Einstein-Maxwell equations. 
Equations (6) then determines the tetrad field only up to local Lorentz transformations 



A(x) 



n ■ 



which retain the condition a M = 0. Hereafter we shall refer to this property of the field 
equations as restricted local Lorentz invariance. 

The structure of the Weintzenbock spaces with spherical symmetry and three unknown 
functions of radial coordinate was given by Robertson [20] in the form 
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(11) 



where the vector 60^ has taken to be imaginary in order to preserve the Lorentz signature for 
the metric, i.e, the functions A and D have to be taken as imaginary. Applying (11) to the 
field equations (Eq.(6)) the author got [19] a set of non linear partial differential equations. 
The solution of these equations has the form [19]: 
First Solution 



If A{R) 



and 

Second Solution 



D(R) = 0, 




B(R) 



1- 



2m 



R= B- 



R 2 ' 



(12) 



// A=l, 



B = 1 



and 



D(r) = 



\J1mr — q' 2 



(13) 



Third Solution 



If A(R) 



RB'Y 



and 



D(R) 



1 



1 - RB' 



n 2 R' 

V + ^( RB ' -2). 
R 4 R ' 




(14) 



It is clear from (14) that the third solution depends on the arbitrary function B, i.e., we 
can generate the pervious solutions of (12) and (13) by choosing the arbitrary function B to 
have the form 



B(R) = 1, and B(R) = J 



1 / 2m q 2 \ , 



(15) 



The associated metric of the three solutions (12) ~ (14) is found to be the same and have 
the form 



dv 

ds 2 = -n(r)dT 2 + —— + r 2 dtt 2 , with n(r) = 
n(r) 



2m q 2 
1 + \ 



(16) 



which is the static Reissner Nordstrom black hole [26, 27]. The form of the vector potential 
A^, the antisymmetric electromagnetic tensor field and the energy- momentum tensor 
are given by [19] 



A(r) 



2v/?r' 



rt 



20Fr 2 ' 



To 



Ti 1 = -T 2 



_T 3 3 



87rr 4 ' 



(17) 



It is of interest to note that the two tensors H^ u and K^ v are vanishing identically for the 
three solutions given by Eq. (12), (13) and (14). It is proved that these tensors are vanishing 
identically for any spherically symmetric solutions [28, 29]. 



3. Singularities 



In teleparallel theories we mean by singularity of space-time [25] the singularity of the 
scalar concomitants of the torsion and curvature tensors. 

Using the definitions of the Riemann-Christoffel curvature tensor, Ricci tensor, Ricci 
scalar, torsion tensor, basic vector, traceless part and the axial vector part [30] we obtain 
for the first solution of (12) 

4g 4 



R 8 



7q 4 - \2RMq 2 + QM 2 R 2 
-2 



R^R,, 



R = 0, 



R 4 (R 2 - 2MR + q 2 ) 



AR A - 12MR 3 + 6R 2 q 2 - AR 3 ^j R 2 — 2MR + q 2 - lORMq 2 



+8R 2 ^jR 2 — 2MR + q 2 M - AR^J R 2 - 2MR + q 2 q 2 + 9M 2 R 2 + 3q 4 



-1 



R\R 2 - 2MR + q 2 ) 



-2R 3 + AMR 2 - 2Rq 2 + 2^R 2 - 2MR + q 2 R 2 



1 2 



-3 ^R 2 - 2MR + q 2 MR + yj ' R 2 - 2MR + q 2 q< 



R 3 - 2MR 2 + Rq 2 - \J R 2 - 2MR + q 2 R 2 



R\R 2 - 2MR + q 2 ) 
+3\jR 2 - 2MR + q 2 MR - 2\] R 2 - 2MR + q 2 q 2 



a^a^ = 0. 



(18) 



The scalars of the Riemann-curvature tensor, Ricci tensor and Ricci scalar of the second 
solution (13) are the same as given by (18). This is a logic results since both solutions 
reproduce the same metric tensor and these scalars mainly depend on the metric tensor. 
The scalars of torsion tensor, basic vector, traceless part and the axial vector part of the 
space-time given by solution (13) are given by 



(2Mr - q 2 )r 4 
-1 

(2Mr - q 2 )r A 



(3g 4 - 10g 2 Mr + 9MV) 
(3Mr - 2q 2 ) 



-1 



(2Mr - q 2 )r 4 



(3Mr - q 



It is clear from (18) and (19) that the scalars of the torsion, basic vector and the traceless 
part of the first two solutions given by Eqs. (12) and (13) are quite different in spite that they 
gave the same associated metric (16)! The singularity of the scalars of Riemann-curvature 
tensor, Ricci tensor and Ricci scalar is given at R — > and this is well know from general 
relativity as we can see from Eq. (18) [1]. As the singularities of the scalars of the torsion, 
basic vector and the traceless part of the first solution (12) are r — > and R 2 — 2Mr + q 2 — > 
the second singularity may has the form R — > M ± \JM 2 — q 2 which is the horizons of the 
static Reissner Nordstrom black hole [31]. 

The singularities of the second solution (19) are given by r — > and q 2 /2M. Now we 
have two solutions reproduce the same metric but the singularity of their space-times are 
not coincide. This is expected of course due to the following facts: 

i) The energy content of these space-times are different [19]. 

ii) The time-space components of the tetrad fields b a , b a ° go to zero as —= at infinity 
[28, 29]. 

iii) Also another interpretations, which may be taken into account as is clear from Eqs. (18) 
and (19) is that the torsion tensor of these solutions is different. As we discussed in the 
introduction that the torsion plays the role of the force, therefore, we may interpret the 
different results of the two torsions given by Eq. (18) and (19) due to the fact that the forces 
of the two solutions are different. 



4. The Stability condition 



In the background of gravitational field the trajectories are represented by the geodesic 
equation 

d 2 x x rx i dx» dx u , . 

!?- + {Mir£- = ' (20) 

dx^ 

where — — is the velocity four vector, s is a parameter varying along the geodesic. It is well 
ds 



know that the perturbation of the geodesic will leads to deviation [1] 

d 2 ( x (x 1 dx^dC ( X ^ dx»dx» n . . 



where ( p is the deviation 4-vector. 

Using Eqs. (20) and (21) in (16) we get for the geodesic equations 



d 2 t 1 . / df 



^0 

ds 2 



= 0, 



d 2 



ds 2 



0, 



(22) 



and for the geodesic deviation 

d 2 (° | j(r)dtdC l 
ds 2 rj(r) ds ds 



0, 



g^C 1 , x „ ,dtd(° . .d(pd( 3 

+ r]{r)r] (r) — 2rr?(r) 



ds 2 
+ 



(is rfs 



cis ds 



rf 2 C 3 2rf0rfC X 



C x = o, 



+ ~p C = o, 



(is 2 \ ds 



where r)(r) is defined by (16), rj'(r) 



ds 2 r ds ds 



0, 



(23) 



plane 

Using (24) in (16) we get 
from (25) and (22) we obtain 



drj(r) 
dr 



and we have consider the circular orbit in the 



9 71 ^ ^ r 
2 ' ds ds 



(24) 



(25) 



r)'(r) 



ds J r(2r)(r) — rr)'(r)) ' 



ds 



2r)(r) — rr]'(r) 



(26) 



The variable s in (23) can be eliminated and we can rewrite it in the form 
d 2 (° r]'{r) dt dC 1 



+ 







d(f> 2 rf(r) dcj) dcj) 

d 2 C , w/ ,dtd(° n . ,d( 3 
__ + jj^ji (r) — — - 2rrj{r) — 

+ 



d(f> d(f> d(j) 
\ (f{r) +v(rW'(r)) [^j - (rj(r) + rrj'(r)) 
d 2 C 2 d( l n 



C x = o, 



d(f> 2 ' d(f> 2 ' r d(p 



+ -- 



(27) 



It is clear from the third equations of (27) that it represent a simple harmonic motion, this 
means that the motion in the plan 6 = n/2 is stable. 



Assuming now the solution of the remaining equations given by 

(° = A ie 1 ^, C 1 = A 2 e^, and ( 3 = A 3 e 1 ^, (28) 
where Ai, A 2 and A 3 are constants to be determined. Inserting (28) in (27) we get 

r 3 m — 6m 2 r 2 + 9mrq 2 — 4g 4 



r 2 (mr — q 2 ) 



> 0, (29) 



which is the condition of the stability for a static spherically symmetric Reissner Nordstrom 
solution. Condition (29) can be rewritten as 

a 2 

r - — > and r - Qm > 0. (30) 

m 



5. Main results 



The main results can be summarized as follows 

1) The singularity problem of the first two solutions (12) and (13) obtained before [19] has 
been studied. The scalars of the torsion tensor, basic vector and the traceless part of these 
solutions are quite different as we can see from Eq. (18) and (19). The scalars have a 
common singularity if r — > 0. Furthermore, the first solution has another singularity if 

r 2 - 2rm + q 2 -> 0, 

while the second solution has another singularity if 



This explains that the structure of the two solutions (12) and (13) are quite different in spite 
that they reproduce the same metric space-time. 

2) The stability condition for the metric of Reissner Nordstrom black hole Eq. (16) is 
derived and is given by Eq. (29). From this condition we can see that. 

i) If r — > the value of (29) is finite. 

ii) If r becomes large then Eq. (30) takes the value r > 6m and r > q 2 /m which is the 
condition of stability for Reissner Nordstrom black hole. 

iii) When q = and if r becomes large then Eq. (32) takes the value r > 6m which the 
is condition of stability for Schwarzschild black hole [30]. The analysis given here for the 
derivation of the stability condition is a straightforward and so simple than that used in the 
literature [32, 33] 
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